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Abstract 

We perform a non-perturbative analysis to the Hamiltonian constraint of the 
lowest-order effective action of the complete Hofava theory, which includes a 
(dilnN) 2 term in the Lagrangian. We cast this constraint as a partial differ- 
ential equation for N and show that the solution exists and is unique under a 
condition of positivity for the metric and its conjugate momentum. We interpret 
this condition as the analog of the positivity of the spatial scalar curvature in 
general relativity. From the analysis we extract several general properties of the 
solution for N: an upper bound on its absolute value and its asymptotic behav- 
ior. In particular, we find that the asymptotic behavior is different to that of 
general relativity, which has consequences on the evolution of the initial data and 
the calculus of variations. Similarly, we proof the existence and uniqueness of the 
solution of the equation for the Lagrange multiplier of the theory. We also find a 
relationship between the expression of the energy and the solution of the Hamil- 
tonian constraint. Using it we prove the positivity of the energy of the effective 
action under consideration. Minkowski spacetime is obtained from Hofava theory 
at minimal energy. 
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1 Introduction 



There are many interesting analyses that can be done in the framework of the Hamil- 
tonian formulation of the Hofava theory pp. Such studies are facilitated by the fact 
that the theory is originally formulated in Arnowitt-Deser-Misner (ADM) variables. In 
this context, the closeness of Dirac's procedure for the preservation in time of the con- 
straints is of central importance for the consistency of the theory. This requires to 
analyze carefully the existence and uniqueness of the solutions of the constraints that 
emerge from Dirac's procedure in the sense of formal problems of partial differential 
equations (PDEs). This is our main goal in this paper. In particular, we shall study the 
existence and uniqueness of the solution of the Hamiltonian constraint of the lowest- 
order effective model of the nonprojectable Hofava theory. The strategy is to regard 
this constraint as a PDE for the lapse function. We shall also study general properties 
of the solutions of the constraints and compare them with known results in general 
relativity. We are interested in analyzing the asymptotic behavior of the field variables 
since, because of the reduced gauge symmetry group, the formulation is different (and 
simpler) to the asymptotic behavior in the canonical formulation of general relativity. 

In the ADM formulation of general relativity, the perturbative Hamiltonian con- 
straint is solved for the scalar mode of the transverse sector of the spatial metric [2J. 
The reduced Hamiltonian is then formulated for the independent modes, the transverse- 
traceless tensorial modes. In turn, the ADM energy is determined precisely by the 
asymptotic behavior of the transverse scalar mode. Therefore, there is a direct relation- 
ship between the Hamiltonian constraint and the energy, the solution of the former is 
the one that contribute to the energy. Since in our analysis we deal with the general 
solution of the Hamiltonian constraint of the Hofava theory, in this paper we also search 
for a relationship between the variable fixed by the Hamiltonian constraint and the en- 
ergy of the Hofava theory. We shall perform the study on nonperturbative grounds. We 
remark that high-order terms in the potential do not contribute to the energy of the 
theory since the energy is given by a surface integral at spatial infinity and any high- 
order term decays faster than 1/r 3 . Therefore, the energy of the lowest-order effective 
theory coincides with the energy of the full theory. 

The Hofava theory can be formulated in two disconnected versions: the projectable 
formulation, in which the lapse function depends only in time and the nonprojectable 
formulation were the lapse function depends both in space and time. Although both pos- 
sibilities are compatible with the gauge symmetry of the theory, the foliation-preserving 
diffeomorphisms, the former cannot be recovered from the latter due to the constraints 
of the nonprojectable case [3]. In the projectable version there is not any version of local 
Hamiltonian constraint (which, in general relativity, is part of the Einstein equations), 
hence general relativity cannot be recovered in the projectable formulation unless radical 
modifications to the theory are imposed, as changing, for example, the gauge symmetry 
group (see [1]). Because of this obstruction we incline ourselves to the nonprojectable 
version, where there is a clear scenario in which general relativity emerges [3]. 

The Hamiltonian formulation of the lowest-order (second-order) truncation of the 
original Hofava theory was achieved by two of us in Ref. [3]. With original we mean 
that the potential is written only in terms of the curvature tensors and cosmological 



2 



constant^, as was formulated originally by P. Hofava in Ref. pQ. However, terms that 
depends explicitly in the spatial derivatives of the lapse function can be included in the 
potential through the acceleration vector a« = ojj In N, which is covariant under foliation- 
preserving diffeomorphisms. These terms were formally proposed for the Hofava theory 
in Ref. Assuming as a principle that a complete, renormalizable, theory must include 
all the terms compatible with the gauge symmetry group, we consider the complete 
Hofava models as those that include the a« terms. 

In Ref. [[3J it was found the rather surprising result that the lowest-order truncation of 
the original theory has a dynamics completely equivalent to general relativity, regardless 
of the value of the constant A that arises in the kinetic term of the model. In a subsequent 
study [6], it was analyzed a model with a R 2 term in the potential, showing the closure 
of the algebra of constraints with the presence of an extra, odd, physical mode. The 
odd dimensionality of the extra mode was previously pointed out in Ref. [7J (see also [8] 
for early evidencejj. 

For the complete theory, the Hamiltonian analysis of the lowest-order effective action 
[T2l[T3"im] showed a consistent structure of constraints. In particular, the model has the 
momentum constraint as the first-class constraint and the Hamiltonian, together with 
the vanishing of the momentum conjugated to the lapse function N, as the second-class 
constraints^]. In addition, the preservation in time of the second-class constraints imposes 
a elliptic PDE for a Lagrange multiplier of the theory. Among the set of constraints, the 
main difference with general relativity is that the Hamiltonian constraint depends on 
the lapse function (as in the case of the original theory with high-order curvature terms 
[6J [7]). This suggest that this constraint can be solved for N on general grounds, since 
it is a linear, elliptic, PDE for N. We note that in (nonperturbative) general relativity 
there is not any plan to solve the Hamiltonian constraint for a specific canonical variable 
that works consistently in general. In this regard the situation looks better in Hofava 
theory, but at the same time elucidating the existence (and uniqueness) of the solutions 
of the Hamiltonian constraint becomes an urgent task for the consistency of the theory. 

We remark that it is only after all the constraints have been solved that one can 
known where are the independent physical degrees of freedom of the theory. This infor- 
mation is crucial in order to study further physical properties of the propagating modes. 
For example, in the perturbative analyses that have been done in the complete theory 
the lapse function is excluded from the independent modes while the extra mode is pa- 
rameterized in terms of the transverse scalar of the metric and its conjugate momentum. 
This procedure is implicitly supported on the assumption that the lapse function can be 
consistently solved from the Hamiltonian constraint, which is the only constraint where 
it gets involved. Any failure in solving the Hamiltonian constraint for the lapse function 
could be the signal of a different independent mode and could change drastically the 
conclusions about the physics of it. Therefore, the rather non trivial structure of the 

4 The model in [3] does not include cosmological constant. 

5 Some perturbative analyses of the original nonprojectable theory can be found in Refs. |9] and 
for the complete theory in [5J [TO] . For a perturbative analysis in the projectable case see, for example, 
Ref. [UJ. 

6 In the complete Hofava theory, unlike general relativity, the lapse function must be considered as 
a canonical variable. 
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PDE for iV that yields the Hamiltonian constraint gives rise relevance to the analysis of 
the existence of its solutions. 

As we mentioned at the beginning, in this paper we elaborate on this point, showing 
rigorously that the solution, in the sense of distributions, exists and is unique. To 
arrive at this conclusion we assume certain condition of positiveness for the gravitational 
variable and its conjugate momentum. This assumption can be seen as the extension 
of the assumption of positiveness of the scalar curvature in general relativity, as we 
are going to see. We shall also study the existence of the solution of the equation for 
the Lagrange multiplier of the theory, in this case the exact solution since the operator 
is simpler. The analysis will also allows us to extract several general properties of the 
solutions for the lapse function and the Lagrange multiplier. We give in advance that the 
evolution equations preserve the leading mode (non-fixed by the boundary conditions) 
of the asymptotic expansion of the initial data; that is, the time derivative of the spatial 
metric has a fall off faster than the metric itself. 

Finally, we shall perform a change of variables by doing a canonical transformation 
that will allow us to reformulate the solution of the Hamiltonian constraint in terms of 
certain combination of fields. This combination depends on N, hence the solution can 
still be used to fix N. The novelty is that the combination is the one that precisely 
yields the energy, much like (perturbative) general relativity. 

The Hamiltonian we will analyze is also related to the famous Perelman energy 
function [13], although the evolution equations are different. It is then relevant its 
analysis also from this different point of view. 

2 Asymptotics in general relativity 

In order to study the similarities and differences between general relativity and the 
Hofava theory in the asymptotic behavior of the field variables, we start with a quick 
summary on the subject in general relativity. We cast the condition of asymptotic 
flatness by demanding that the coordinate system becomes Cartesian at spatial infinity. 
This approach is convenient for us because of the connection with the problem of solving 
the PDEs that arise in the canonical formulation of the Hofava theory. 

In the standard (ADM) canonical formulation of general relativity one starts by 
splitting the space-time into spacelike submanifolds and the time direction. To this end 
one needs a concrete notion of time, hence a global time function t is introduced. On 
the basis of this function the space-time is regarded as a foliated manifold M = E x K, 
where E represents the three-dimensional spacelike slices of constant t and one also 
assumes that it is endowed with a Riemannian metric gry. Since in this work we are 
interested in gravitation as a local theory (specially for isolated systems), we consider 
£ to be a complete, connected, asymptotically flat Riemannian manifold. 

Let n M be the future-directed, unit vector field orthogonal to £ with respect to the 
space-time metric. Being n M orthogonal to E, its action over t does not yield zero. In 
addition, the future-directed condition restrict n^d^t to be positive. The lapse function 
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N is determined by the foliation and the space-time metric by the formula 



iV = — - — . (2.1) 

Let V 1 be a vector field satisfying 

t^d^t = 1 . (2.2) 

One may decompose t M into its normal and tangent components with respect to E. By 
using (12. ip in condition (12.21) . one obtains that t M takes the form 

tP = Nn» + N» , (2.3) 

where N 11 , the shift function, is the part of i M lying over S. If one chooses t as the time 
coordinate and also takes t M such that the time coordinate is adapted to it, t^d^ = d t , 
then the space-time metric takes the well-known ADM form, 

ds 2 = (-N 2 + NiN^dt 2 + 2N i dx i dt + g ij dx i dx j , (2.4) 

where Aj = gijN 1 . 

It is important to define two more objects associated to the foliation. One is the 
extrinsic curvature K^, which basically measures the deviation of n M from itself under 
parallel transport along S. Since P^ a = 5^ a + n p n a is the projector to S, K^ v = 
Pn a Pv^ a np. The spatial components of K^ v constitute a tensor over S and are given 
by 

^' = ^(^-V(^))- (2-5) 
They essentially play the role of conjugate momentum of in the canonical formulation, 

^ — — g ij K) . (2.6) 

The other object is the acceleration vector, 

a i = d i laN, (2.7) 

which is the gradient of along E. This vector plays a key role in the nonprojectable 
formulation of the Hof ava theory since it enters in the potential [5] . 

The asymptotic behavior of g^ and can be conveniently established in terms 
of the coordinate system. For r — > oo, where r = vW and x % become Cartesian at 
spatial infinity, an asymptotically flat metric tends to the Euclidean metric in the form 

9ij = 5 tJ + 0(r~ l ) , (2.8) 

whereas the extrinsic-curvature tensor falls off at the speed 

K l3 = 0(r- 2 ) . (2.9) 

Actually, these are consequences of the formal definition of asymptotic flatness, which 
roughly consists of the requisite that the metric can be completed at infinity in the same 
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way that Minkowski space-time can be completed by a conformal infinity (see, e. g., [16] 
and see [17] for a definition of asymptotic flatness for S that involves condition (12. 8p 
directly) . 

Condition (12. 8p has important consequences in the way energy is defined in general 
relativity. In this theory there is not a universal expression for a local energy density, 
such that the energy could be formulated as a volume integral. Instead, the energy, as 
well as the momentum, are formulated as surface integrals at spatial infinity. Specifically, 
for asymptotically flat space-time, the ADM energy and momentum of the gravitational 
field are respectively defined by 



■^ADM 



j dyijlj),!!, : - dig j:j ) , (2.10) 
P1 DM = -2<taZ^. (2.11) 



It turns out that £ ADM has physically acceptable values, that is, finite and different from 
zero in general, if 9ij behaves at spatial infinity as indicated in flZEjfl. Similarly, the 
decay (EHD yields 

tt« = 0{r- 2 ) , (2.12) 

such that the configurations have finite and in general different from zero momentum. 

It is worth stressing that the ADM energy (I2.10p of all configurations is determined 
only by the term of order r _1 once gij is expanded at the asymptotic limit (e. g., the 
mass parameter in Schwarzschild space-time). Therefore, when performing calculus of 
variations, variations of order bg^ = (9(r _1 ) allow to compare between points of the 
phase space of (generically) different energy. On the contrary, any variation of faster 
decay preserve the value of the energy. 

The Hamiltonian of general relativity is given by 

H GK = J d 3 x(N'H GK + N{H GR ) + E Ar>M , (2.13) 

U GK = -^(tt^- - Itt 2 ) - y/gR, (2.14) 

W Gn ee -2V^\ (2.15) 

The ADM energy is included in the Hamiltonian in order to get the correct evolution 
equations under arbitrary variations, which are considered to be bg^ = (9(r _1 ) asymp- 
totically [18] . For variations with this asymptotic fall off, the variation of the curvature 
term, 5(—y/g~R), gives rise to a nonvanishing boundary contribution that is exactly 
—SE AVM , hence it cancels completely the variation of the boundary term in (12. 13ft . If 
the -Eadm term was not included in the Hamiltonian, then the boundary contribution of 
5(—^/g~R) would be an obstruction to get any evolution equations. This argument is an 
elegant way to justify the ADM energy of the gravitational field [IB] . Since the bulk part 
of the Hamiltonian is a sum of constraints, H GR = ^ 8 GR = 0, when the Hamiltonian is 



7 For this statement we assume that the spatial derivatives of any mode of order r P, (3 > 0, are of 
order r^^ +1 \ 
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evaluated on any physical configuration it yields the corresponding value of E ADM . More- 
over, the E AVM term of the Hamiltonian can be obtained from the very covariant action 
by requiring its differentiability under arbitrary variations, Sg^ = (9(r _1 ) asymptoti- 
cally [12]. For variations with faster decay there are not boundary contributions coming 
from the bulk part of the Hamiltonian and SE ADM is automatically zero. 

The canonical formalism is particularly suitable to formulate initial-data problems 
in general relativity. As can be seen from (I2.13p . in general relativity N and iVj can be 
dropped out from the phase space since they are the Lagrange multipliers of the con- 
straints % and TV. These are first-class constraints, hence N and Ni can be determined 
by gauge fixing the symmetries of general relativity. Therefore, the initial data consists 
of an asymptotically flat set (£, , 7T'- 7 ) given at the initial time and subject to the 
constraints (I2.14p and (I2.15p . The evolution equations 

2N 

gij = —(7r ij -y ij 7r)+2V ii N j) . (2.16) 

9 AT 

jc ij = -N^g (RV - y*R) - — [tt V - \^ - ^'(Afci - l^r 2 )] 

+ v /#(V i <W - g ij V 2 N) - 2V fe iV ( V )fc + V k {N k 7T ij ) , (2.17) 

determine the evolution in time of the initial data once the lapse and shift functions 
have been chosen; that is, the evolution is given in a particular gauge. The evolution 
equations automatically preserve the constraints (12.141) and (12.151) since these are first- 
class constraints. The initial data (£, g^, k^) contains twelve functions over E; the four 
constraints and the four components of an arbitrary diffeomorphism reduce the number 
of independent initial functions to four. 

The standard assumption [18J for the asymptotic behavior of the lapse and shift 
functions is that, in the coordinate system that becomes Cartesian at infinity, they 
behave as 

N = l + O^ 1 ), (2.18) 

N t = 0{r- 1 ). (2.19) 

These conditions are read from the Schwarzschild solution and, together with ( 12.81 ) 
and (I2.12p . allow for the differentiability of the Hamiltonian. However, the structures 
(l2.18H2.19l) are not the most general ones since the phase space must include those con- 
figurations obtained by the action of the symmetry group of the theory. Coordinate 
transformations that do not preserve the foliation do not preserve (12.1 8H2.19]) . In par- 
ticular, there is the possibility of performing a coordinate transformation that involves 
the action of the Poincare group at spatial infinity. For the case of a pure, infinitesimal, 
boost at spatial infinity the lapse function takes the form 

N = fix 1 + Oir- 1 ) , (2.20) 

such that it diverges linearly. This is important in order to establish the compatibility 
of the asymptotic behavior of the field variables with the evolution equations. From 
(12.161) and (12.171) and setting iV, = one has, asymptotically, 

g l3 = 0^) , tt« = 0(r- 2 ) . (2.21) 
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Therefore, the decays (12.81) and (I2.12p are preserved. 



3 Hamiltonian and constraints of the Hof ava theory 

In this section we give a summary of the Hamiltonian formulation for the lowest-order 
effective action of the complete nonprojectable Hof ava theory [121 1131 El- The gauge 
symmetries of the Hofava theory, which are the foliation-preserving diffeomorphisms, 
correspond to coordinate transformations of the kind 

t' = t'(t), x' = x'(t,x), (3.1) 

This class of transformations is a subset of the ones allowed in general relativity. The 
theory is written in terms of the Arnowitt-Deser-Misner variables, with the novelty that 
the Lagrangian of the complete theory depends on the spatial vector a, defined in (12.71) 

B- 

The second-order action (without cosmological constant) of the complete theory is 
given by 

S = J dbd z x s /gN{G ijkl K ij K kl + R + aa^) , (3.2) 
where a is an arbitrary coupling constant and 

K* = 2^(fo - 2V (< JV0) , (3.3) 

Whenever A 7^ 1/3, the inverse of G l ^ kl is given by 

Gijki = -(dikdji + 9u9jk) + _ ^^ 9ij9ki ■ (3.5) 

The corresponding Hamiltonian takes the form of a sum of constraints plus boundary 
terms [HI E31 E] , 

H = J d 3 x(Nn + NiH* + 0-0) + E ADM - 2a$ N , (3.6) 

H = G lj ki 7 ^^ + ^g{-R + a{2V t a l + a t a l )) ) (3.7) 

W ee -2Vj7i ji + 0N , (3.8) 
<5> N ee idZidiN, (3.9) 



where <fi is the momentum conjugate to N and a is a Lagrange multiplier. The momen- 
tum constraint W = is the first-class constraint of the theory whereas = and the 
Hamiltonian constraint % = are the second-class ones. The last term in the momen- 
tum constraint (13. Sp is included in order to this constraint becomes the full generator of 
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spatial diffeomorphisms, since the combination <pdiN serves as generator of these diffeo- 
morphisms on N and (p. Notice that the dependence of the momentum constraint on N 
vanishes on the constrained phase space, hence the Hamiltonian constraint is the only 
one where N gets involved once all constraints are imposed. As in general relativity, the 
ADM energy is included as the boundary term needed to obtain the equations of motion 
under variations bg^ of order (9(r _1 ) asymptotically. Similarly, the flux of N at spatial 
infinity, $/v, cancels a nonzero contribution coming from 5(2aN'Via l ) for variations of 
iV with asymptotic decay 5N = (9(r _1 ). 

To preserve in time the second-class constraints a second-order, linear, PDE for a 
must be imposed. This equation was found in Ref. [H] in the form 

V 2 (a/iV) + a% (cr/N) = J, (3.10) 

J = -T^Vv^VV) + ^-G t3kl [V\a^ kl ) + 2aVvr fc '] . (3.11) 
y9 \J9 

where 7 = (1 — A)/(l — 3A)a. We can be brought this equation to the form 

diiN^d^a/N)) = di (- 7 iV 2 VV + N 2 g l klm a k n lm ) , (3.12) 

such that the source term is also an exact divergence. 

In Hofava theory the initial data is less restricted than in general relativity, but 
the propagation of the constraints requires the knowledge of the master solution of 
Eq. (I3.10p valid at each instant of time (see Ref. [13]). The shift functions Ni are still 
the Lagrange multipliers of the momentum constraint W, hence they are again fixed 
by a procedure of gauge-symmetry fixing. The general plan, for which we are going to 
study its consistency in detail in the next section, is to cast the Hamiltonian constraint 
% = as a PDE for the lapse function N. Therefore, the initial data consists of the set 
(E, gij, tt^) subject only to the momentum constraint H 1 = 0. At the initial time, the 
lapse function iV must be solved from the equation TL — 0. Then, the evolution of the 
initial data is dictated by the canonical equations of motion of the theory, 

kl 

g l3 = 2NQijki—-= + 2V(jiVj) , (3.13) 
\/9 



+^/gN [VV + (1 - a)aV - g ij (V k a k + (1 - \a)a k a k )] 

-2V k N (i 7i j)k + V k (N k 7i i3 ), (3.14) 
A> = (T + rfdiN, (3.15) 

once iVj has been gauge fixed! As we have mentioned, for the whole propagation 
it is assumed that Eq. f)3.10p is systematically solved for a at each instant of time, 
starting with the initial time where the initial data is given and N(0) is solved form the 



is automatically zero because of the constraints of the theory. 
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Hamiltonian constraint "H(O) = 0. This requisite is a consequence of the second-class 
constraints of the theory, unlike general relativity. In the next section we also are going 
to study the existence and uniqueness of the solution of Eq. (13.101) . Now the initial 
functions over E are restricted by the three components of the momentum constraint 
and there are three independent components of a spatial diffeomorphism. This leaves 
us with six independent initial functions over E. 

4 Existence and uniqueness of the solution of the 
Hamiltonian constraint 

We will now discuss the existence and uniqueness of the solution of the Hamiltonian 
constraint (13. 7\i and afterwards of the solution of Eq. (I3.12p . We shall consider global, 
nonperturbative, arguments which will be also useful in the constrained problem in the 
full Hofava theory. First, it is convenient to rewrite the constraint (13.71) in the following 
way. Let iV be continuous with right and left derivatives, then 

N = sgn(N)\N\ . (4.1) 

It yields 

diN = sgn(N)di\N\ (4.2) 

because <9jSgn(iV) is a Dirac delta distribution on the direction orthogonal to the sub- 
manifold N = 0, but is is multiplied by \N\ which is zero on the support of the delta 
distribution. We then have 

m o t \N\ 



N \N\ 



diln\N\. (4.3) 



The equation ti — 0, where H is given in (13. 71) . may be formulated as an equation for 
\N\, 

-\7 2 ^/\N\ + ^(R- Gm^-j VW\ = o. (4.4) 

Hence, besides proving the existence and uniqueness of the solution one has to show 
positiveness of it. It will be convenient to write v/pV] = 1 + u, then Eq. (I4.4p becomes 



i / ijkl\ 

- d^g^dju) + Gu = -G, G = — (ygR - & jk i^=-) ■ (4.5) 

Let u be a scalar function on E. We introduce the space C' 1 (E) according to 

C^E) = { u e C^E) : u -»■ when r oo} . (4.6) 

C 1 (E) is the space of functions that have continuous derivatives. It is invariant under 
diffeomorphisms on the slice E. 

We denote, for u and v in C 1 (E), 

{u, v)= [ Sxi^ggVdiudjV + Guv) , (4.7) 
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where G is defined in (I4.5p . The map (14. 7p is a bilinear functional on E. In the definition 
of this map we assume that R is integrable over E. If G > on E, then (u, u) > and 
(u,u) = if and only if u = 0. In fact, (u,u) = implies J d 3 x^/gg 13 diudjU = 0, hence 
u = constant on E and the boundary condition implies that u = on E. Consequently, 
if G > on E, (u,v) defines an internal product on (7 1 (E). 

Let "H 1 denote the Hilbert space obtained by the completion of C ,1 (E) with respect 
to the norm induced by the internal product (14 .7p . We can now prove the following 
Proposition 1: Given a complete, asymptotically flat Riemannian manifold E, mo- 
menta 7r y over it and values of the parameters A and a satisfying G > 0, there always 
exists in "H 1 a unique weak solution u of Eq. (14. 5ft . 

Proof : In order to prove this statement may use the Lax-Milgram theorem. However, 
in this case it is a straightforward application of the Riesz theorem. In fact, for any 
v G Ti 1 , jj^Gv is a linear, continuous, functional with domain "H 1 . We notice that Gv 
is integrable at infinity. Riesz's theorem ensures the existence of a unique w G "H 1 such 
that 

(w,v) = -J Gv (4.8) 
s 

for all v G Ti 1 . But Eq. (14 .8p is exactly Eq. ( I4.5P for w when the latter is applied as a 
distribution to a test function v G 'H 1 with compact support. In fact, ^/gg^diwdjV and 
Gv are scalar densities on E that are integrable at infinity for all v G H 1 , in particular, 
for C°° test functions of compact support. We then have a unique solution of Eq. (14.51) 
in the sense of distributions. 

Remark: The assumption G > is essential to prove existence and uniqueness of the 
solution. If G is bounded from below, G > —C, C > 0, then in a generic situation 
— V 2 + G has a nontrivial kernel K. In this case Eq. (14. 5p has a solution if and only if 
G is orthogonal to the kernel K. The solution, if exists, is the addition of a particular 
solution plus the zero modes. It is then not unique and it may even not exists. 

We can improve the result, as usual in the theory of elliptic operators, to show that 
the solution is C°° if gij and tt^ are C°°. To do so, we first characterize the asymptotic 
behavior of the solution. 

Proposition 2: Given an asymptotically flat Riemannian manifold with C°° metric g^ 
and momenta 7r y- , then the asymptotic solution of Eq. (14. 5 j) exists, it is C°° for large 
enough r and its asymptotic behavior is u = 0(r~ l ). 

Proof : The asymptotic version of Eq. (14.51) is equivalent to perturbate it around the 
Minkowski solution. Since the momentum 7r* J enters quadratically in Eq. (14. 5ft . it diss- 
apears from this equation at linear order. Hence, for our purposes it is enough to deal 
with the perturbative metric and u. Let us expand the metric at the asymptotic limit, 

9ij = 8%j + hij , (4.9) 

such that hy = (9(r _1 ). At linear order in h^, the scalar curvature takes the form 

R = d ij h ij -d ii h, (4.10) 

where h = hjj. We use the shorthand 9y = didj, and so on. As usual in perturbation 
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theory, we decompose the perturbative metric into transverse and longitudinal parts, 

hij = hf? + \\5ii - d^dkuY 1 }^ + dihf + fyti + d^duuY 1 ^ , (4.11) 

where hj^ = 0, dihj^ = 0, dihf = and the boundary conditions are used in the 
definition of (dkk) 1 - All of these variables are determined uniquely in terms of h%j. The 
scalar of the longitudinal sector is given by 

h. L = (d kk )- 1 d ii h ij , (4.12) 

and starting from this the other variables are solved: the longitudinal vector hf is 
solved from the divergence of Eq. (14. lip , djhij = djjhf + dih L , the transverse scalar 
is obtained from h = h T + h L , and hf^ can be solved directly from (14. lip . After 
decomposition (14. lip , the linearized scalar curvature (I4.10p is determined only by scalar 
of the transverse sector, 

R=-diih T . (4.13) 

We now expand the Eq. (14.51) at r — > oo. We notice that among the terms that depend 
on u the dominant one is the Laplacian. In the right hand side, the dominant term in 
G is R. Thus, for large enough r, the Eq. (14.51) takes the linearized form 

d ii u = d ii h T . (4.14) 

The solutions to this equation exist and are given by u = h T plus an element of the 
kernel of da. The boundary condition u\oo = fixes the solution u = h T . Clearly, the 
solution is C°° provided that the metric is C°° and near spatial infinity it behaves as 
u = h T = Oir- 1 ). 

Now we are going to show that, as a consequence of Proposition 2, the weak solution 
of Proposition 1 is asymptotically C°°. 

Proposition 3: Given a complete, asymptotically flat Riemannian manifold E with C°° 
metric and conjugated momenta such that G > on E, the weak solution of Proposition 
1 is C°° on E and has the asymptotic behavior of Proposition 2. 

Proof : For the proof we use the elliptic regularity theorem. Under the assumptions 
G G Hl° c for all s, hence the weak solution belongs to H l °+ 2 f° r a ^ s - We then conclude 
that u G C°° and it has the asymptotic behavior of Proposition 2. 
Proposition 4: Under the same assumptions of Proposition 3 and assuming also that 
E is path connected, the solution of Eq. (14.51) satisfies 1 + u > 0. 

Proof : We proceed by contradiction. Let us suppose that at some point p on E the 
solution satisfies u(p) < —1. At spatial infinity u — > 0, hence, by continuity, on a path 
form p to infinity there exists a point at which u = —1. We denote w -1 (a) the set of 
points of E at which the solution is u = a, and also M a = {q G E : l+u(q) < a}. If —1 is 
not a critical value of u, that is, | \diu\ \ ^ for all q G u~ x {— 1), then the implicit function 
theorem ensures that n _1 (— 1) is a submanifold of E. Moreover, m _1 (— 1) is closed and 
bounded. Since E is complete then M is compact with boundary dM = u~ l (— 1). 
From Eq. (14. 5 j) we obtain on M : 

J dariig'^il + ^djil + u) = J d 3 x [y/gg ij di(l + u)dj(l + u) + G(l + u) 2 ] . (4.15) 

dM M 
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But at 8Mq we have 1 + u = 0, thus 

J d 3 x ygg ij di{\ + u)dj{\ + u) + G(l + uf] = , (4.16) 

Mo 

which yields 1 + u — constant. Since at 8Mq it holds that 1 + u — 0, then 1 + u = 
everywhere on M . If —1 is a critical value of u, we consider m _1 (— 1 + e) and use the 
continuity of the integral in M e as e — > to conclude, by the same arguments, that 
1 + u = on M . We thus reach a contradiction with the assumed condition u{jp) < — 1 
in each case. Therefore, we have 1 + u > and we can identify \N\ = 1 + u. 
Proposition 5: Under the assumptions of proposition 4, the solution of the Eq. ( 14.51) 
satisfies |JV| < 1. 

Proof : Let us suppose that \N\ = a > 1 at some point p on E. Then the set M a = 
{p G £ : |iV| > a} is closed and bounded and, since £ is complete, it is compact. Then 
Eq. f !4.4p implies 

I J mg ii V90j\N\ = j(V99 lJ d^\d^\+G\N\). (4.17) 

dM a M a 

That is, the outgoing flux of d{\N\ through dM a is positive. However, from the definition 
of M a we must have that if a < b =^ C M a , which contradicts the flux relation. 
Consequently \N\ < 1. 

Remark: The solution for iV may have +1 or —1 as asymptote. If £ has several 
disconnected sectors at infinity we may impose the +1 or —1 value as the asymptotic 
condition on each of the sectors. The situation is then very different from what occurs in 
the large- distance effective action studied in Ref . [3] , the XR- model, where the constraint 
for iV is 

- diiVg^djN) + ^RN = . (4.18) 
In this case, the conservation of the Hamiltonian constraint, 

7[ ij 7l kl 

Qijki^^- VdR = 0, (4-19) 
V 9 

implies it = 0, which, when combined with the Hamiltonian constraint, implies that the 
scalar curvature is nonnegative, 

R = ^L>0. (4.20) 
9 

All the propositions we have considered here apply also for this case. The assumption 
G > is now replaced by R > 0, which now is a consequence of the constraints of the 
theory. The solution for iV is then positive and satisfies < iV < 1. Note that in general 
relativity, when the tt = gauge is not imposed, the condition R > is not implied by 
the constraints in general. 

Remark: If we take the asymptotic boundary condition |7V| — > at spatial infinity 
then \N\ = on the whole £. Hence this boundary condition has to be dismissed. We 
may take \N\ — > constant 7^ and everything follows in the same way as for \N\ — > 1. 
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Remark: The net outgoing flux of di\N\ is 

\ f n^^d 3 \N\ = J (^di^/lNldjVWl + GlNl) >C. (4.21) 
OS s 

If we take iV — > 1 at infinity then N = \N\ and the flux of diN is always positive. We 
notice that ypVj is not in i-L 1 , hence the right hand side of Eq. (I4.2ip is not a norm. 

The minimal flux condition occurs for N = \N\ = 1. In this case the constraint (|4.4p 
reduces to the Hamiltonian constraint of the XR model (j4.19j) . that is, G — 0. Since the 
evolution equation for N is N = a, Eq. A3. 1 2|) reduces to the corresponding equation 
7i = of the XR theory. From the evolution equation for tt we get R = 0. From the 
Hamiltonian constraint (I4.20p we then get tt^ = 0. We then get Rij = and the three 
dimensional slices £ are flat. Since N = 1 we then regain the Minkowski space-time at 
the minimum flux condition. 

Proposition 6: Under he assumptions of proposition 2, the exact solution of Eq. (13.101) 
exists and is unique. It is C°° and its asymptotic behavior is C(r~ 2 ). 
Proof : If N ^ 0, Eq. ( 13. 10ft expresses an elliptic operator acting on a/N equal to a 
source term. The operator di{N yjgg^ dj ■) is strongly elliptic on the complement of the 
iV = set, which we denote by M . M is compact, on its complement M the solution 
a/N of the Eq. ( 13.101) with boundary condition a/N = on <9£ exists and is unique. 
At the asymptotic limit r — > oo we expand perturbatively the field variables as in ( 14. 9ft . 
also expanding the momenta around zero, 

7r ij = p ij . (4.22) 

At linear order Eq. (I3.12p takes the form 

duo- = -jdup , (4.23) 

where p = p n . With the boundary condition, the unique solution to this equation is 
a = — 7p, which is C°° and of order (D(r~ 2 ). 

Now we are in position to establish the differences in the asymptotic behavior of the 
field variables between general relativity and Hofava theory. The main point is that in 
the nonprojectable version of the Hofava theory the lapse function is not a Lagrange 
multiplier that could be fixed by appealing to any gauge symmetry. Instead, iV is 
determined by the Hamiltonian constraint (14.41) . Therefore, the asymptotic behavior of 
N is dictated by the rest of canonical variables. We have seen that, near spatial infinity, 
the Hamiltonian constraint (I4.4p implies 

N = 1 + Oir' 1 ) . (4.24) 

This decay cannot be changed by any coordinate transformation between the kind al- 
lowed in Hofava theory ( 13. ip . This is in contrast with general relativity, where the 
asymptotic behavior ( 12. 18ft of iV is understood to hold in a specific coordinate system 
that becomes Cartesian at infinity, but it is changed by general Lorentz transformations 
at infinity, in particular by boosts. In Hofava theory the decay (14.241) holds for any 
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asymptotically flat configuration and it is preserved by translations and spatial rota- 
tions at infinity, which are the only symmetries of the Minkowski space allowed by the 
theory. 

As a consequence, the evolution equations (13. 13|) and ( I3.14|) dictate that the time 
derivatives of the metric and momentum decay as 

9ij = 0(r- 2 ) , tt« = 0(r" 3 ) . (4.25) 

These are faster than the decays of g^ and tt 1 ^ themselves at spatial infinity; in contrast, 
again, with the observed in general relativity. This feature also holds for N, which can 
be seen from combining the decay of a, a = 0(r~ 2 ), with Eq. ( 13. 15j) . yielding 

A> = 0(r- 2 ) (4.26) 

asymptotically. These results imply that in Hofava theory a part of the initial data 
is always preserved by the evolution equations: the asymptotic modes of order r _1 for 
<7ij(0) and r~ 2 for 7r y '(0). Similarly, the mode r _1 of N(0) is preserved. 

Finally, the results ( I4.25P and (I4.26P forces us to reconsider the functional variations 
admissible by the Hofava theory, since these variations must belong to the subspace of 
the fields themselves. The most general asymptotic behavior for variations that preserve 
(14725]) and (jPBj) are 

8 9ij = ©(r" 1 ), <5tt« = 0(r- 2 ), 5N = 0{r~ x ) , 

(4.27) 

5 9tJ = 0(r- 2 ), 8ttV = C(r" 3 ), 5N = 0{r- 2 ) . 

5 The energy of the theory and the solution of the 
Hamiltonian constraint 

We consider the conformal transformation 

9%i = e^gij, (5.1) 
n ij = e -Vfti3 f (5.2) 

where the Riemannian metric g^ is subject to 

det^ = l (5.3) 

(this can be imposed in any coordinate system). We also make the covariant decompo- 
sition of 7T 1 - 7 into its traceless and trace parts with respect to the metric cjij, 

* ij = 4 j + , (5-4) 

where g 1 ^ is the inverse of g^, n = ^7r u and 7r^ is traceless, g%p: % ^ = 0. These variables 
will allows us to manage the solution of the Hamiltonian constraint of the theory ( 13.71) 
and also to connect the solution with the expression of the energy, as we are going to see. 
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The transformation (j5.HI5.4p is canonical in the sense that the kinetic terms preserve 
their canonical structure. Indeed, by imposing condition ( 15. 3p . it is straightforward to 
obtain 

d 3 x ir^gij = J d 3 x ("K^g^ + fcip) . (5.5) 

E 

Let us start with the computation of the energy. It is given by the value of the 
Hamiltonian ( 13. 6 p on each configuration. Since the bulk part of the Hamiltonian is a 
sum of constraints, it does not contribute to the energy of any physical configuration. 
Thus, the energy is given by the values of the boundary terms, 

E = E ADM - 2a$ N . (5.6) 

This feature will hold even when high-order terms will be included in the potential of 
the theory. Moreover, since any high-order term has a decay faster that 1/r 3 at spatial 
infinity, they do not contribute to the energy of the theory. Therefore, expression ( 15 .6p is 
actually the energy of the full theory. This one of the aspects that give relevance to the 
lowest-order effective theory. In Ref. [20J it was formulated a positive-energy theorem 
for the complete nonprojectable Hofava theory. One of the hypothesis of the theorem 
is to impose the condition tt = 0. 

After the transformation (l5.1N5.4p . the ADM energy takes the form 



E ADM = j) dE^ {djcpc/ij + djjjij - dicpgjj + d^jj) 



(5.7) 



We may further refine this expression if we assume the asymptotic behavior g^ = 5ij + 
0(1/ r) and <p = 0(1/ r). It follows that at spatial infinity g~jj = 3 and also, by combining 
with condition ( 15. 3ft . di§jj = (det g)~ l di det g = 0. We may also fix a partial gauge in 
which the longitudinal part of vanishes, 

0^ = 0. (5.8) 

To preserve in time this partial gauge-fixing condition one must combine it with the 
evolution equation (13.131) . The resulting equation can be interpreted as the condition 
fixing iVj. After these considerations the ADM energy (15. 7p takes the form 

E ADM = -2 j> d^idi^p , (5.9) 

which only depends on (p. Finally, we get that the energy of the Hofava theory is given 
by the surface term at spatial infinity 

E = —2 j) dY^di (alnN + Lp) , (5.10) 

where for convenience we have written diN as di \nN at spatial infinity. Thus, we obtain 
that the energy of the Hofava theory is determined by the variables iV and ip in terms 
of the combination £ = a In N + (p. 
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Now let us study the Hamiltonian constraint (13. 7\i . By performing the change of 
variables (j5.lll5.4p into this constraint we obtain 



U = e- 3 ^ 2 (np TtJ + (3 - 9A)- 1 7T 2 ) 

+ e^ /2 (-R + 2V 2 £ + acua* + aun&<p + \dnp& 



(5.11) 



In this expression indices have been raised/lowered with g^ and R, V refer to cjij. From 
the 7r 2 term of this constraint we notice that, in order to the mode (p, it) has a good 
physical propagation, the value of A must be restricted by A < 1/3, which is one of the 
admissible ranges found in Ref. [5]. 

As we have seen, the combination £ = a In iV + (p is the variable that determines 
the energy of the Hofava theory. On the other hand, in Hof ava theory the Hamiltonian 
constraint is an equation for N. These facts motivate us to solve the Hamiltonian 
constraint for N in terms of the composed variable £. This can be achieved since in this 
constraint we can combine all the terms that depend on N in a single elliptic operator 
acting on where / is proportional to £. Indeed, it is straightforward to check the 
identity 



diie^gVdje^ 2 ") = — e ^e i/2a 



1 

4a 



2V 2 £ + aa^ + aat&ip + + , (5-12) 



where = (a — 2) /2a. With this identity the constraint H = can be formulated as 
the following equation for £: 



-di(e^g ij d^ /2a ) + GS' 2a = . 



G 



4a 



e~ 2 * (7ttj>7CTij + (3 - 9A)- 1 7T 2 ) - R - pdtf&tp 



(5.13) 
(5.14) 



Note that G depends on the canonical variables §ij, it^, V 9 an d tt> but does not depend on 
N and does not have explicit dependence on £. If \ip\ is bounded, the composed metric 
e^g 1 ! arising in the differential operator —di(e^ v g l Wj •) is positive definite. Thus, the 
operator is coercitive. Under the assumption G > 0, where G is given in (I5.14p . all the 
results of Sec. 4 about the existence of the weak solution hold. Therefore, we may solve 
Eq. (I5.13P for iV if we solve it for the combination £, which is the only place iV where 



arises. 



The Hofava theory may be formulated in terms of g~ij, 7r* 3 ', N. The relation between 
the Hamiltonian constraint formulated in terms of £ and the expression of the energy in 
terms of the flux of yields directly the positivity of the energy of the Hofava theory. 
In fact, from (15.131) we have 



E 



2 <h dT^idii = -4a / d 3 x 



1 



4a 2 



e t/« e P<P§iJd i £d£ + Ge 



(5.15) 



and the right hand member is positive assuming a < 0. We remark that the condition 
G > is imposed to ensure the existence of of the solution for £, and, as a consequence, 
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it leads to the positivity of the energy. Moreover, we have shown in a previous remark 
that N = 1 implies, using the field equations, that gij = Sij and tt^ = 0. Since N = 1 
we end up with the Minkowski spacetime, hence <p = and £ = 0. Consequently, N = 1 
is a configuration of minimal energy, E = 0. It remains open the question if there are 
other configurations of minimal energy. 

The condition a < implies the same sign in the curvature term R and the a^a 1 
term in the Hamiltonian (13. 6p . These terms, together with the factor N, correspond 
exactly to the Perelman energy function 

[(R + VJVifie-fdV, (5.16) 
s 

if we take / = — IniV. The term Qijki 7r '^' i n A3.7[) implies a different dynamics from 
the Ricci flow in Perelman's work. 

6 Conclusions 

In this paper we have shown rigorously the existence and uniqueness of the solution 
of the Hamiltonian constraint of the lowest-order effective action of the complete non- 
projectable Hofava theory. The constraint is solved for the lapse function N and it is 
subject to Dirichlet boundary conditions at spatial infinity. For the proof we have as- 
sumed a kind of positivity condition on the spatial metric and its conjugate momentum. 
We have shown that the Hamiltonian constraint is a PDE for the absolute value of the 
lapse function and, under the same positivity condition, we have obtained a condition 
of consistency: the equation preserves the non negativity of the absolute value. We 
have also found that the lapse function is bounded from above, reaching its maximum 
value at spatial infinity. Of course, it would be interesting to study the existence in the 
general case, without imposing the positivity condition on the spatial metric and its 
conjugate momentum. 

The fact that in the complete Hofava theory, unlike general relativity, the lapse 
function is determined by other field variables allowed us to determine its asymptotic 
behavior at spatial infinity, finding the behavior N = 1 + (9(r _1 ) in any of the allowed 
frames of the theory. As a consequence, the evolution equations yield g = 0(r~ 2 ), a 
(unexpected) behavior different to general relativity. In particular, this implies that the 
evolution equations preserve the mode of order r -1 of the initial data g^it = 0). In 
addition, the variational calculus in the complete Hofava theory must be defined such 
that g = 0{r~ 2 ) is preserved. 

A nice feature of the theory we have found is that, after a conformal transformation, 
the Hamiltonian constraint can be solved in terms of the variable that exactly gives the 
energy of the theory. This is similar to the perturbative ADM Hamiltonian analysis 
of general relativity, where the Hamiltonian constraint is solved for the scalar mode 
of the transverse sector and the energy is defined by the asymptotic behavior of this 
mode. Using this conformal transformation we proved the positivity of the energy of 
the large- distance Hofava theory, under the assumption G > 0. The minimal energy is 
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obtained when N — 1, in this case the theory describes Minkowski spacetime. Although 
the Hofava theory is not invariant under Lorentz tansformations, Minkowski spacetime 
is recovered under the minimal energy condition. 
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